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NOMENCLATURE 

A  a  constant  of  integration 

b  dimensionless  boxindary  layer  thickness  (6f/D) 

B  a  constant  of  integration 

C  solute  concentration  in  liquid  at  distance  x  from  interface,  moles/volume 

C  solute  concentration  in  liquid  at  eqtiilibriim  with  solid  phase, 
a 

moles/volume 

C  eutectic  composition,  moles/volume 

d 

C,  bulk  liquid  concentration,  moles/volume 

C  initial  solute  concentration  in  liquid,  moles/volume 
o 

C  solute  concentration  in  solid  at  equilibrium,  moles/volume 
s 

D  coefficient  of  diffusion 

e  base  of  natural  logarithm 

f  rate  of  freezing  (distance/time) 

g  fraction  frozen 

g  fraction  frozen  of  pure  component 

G  temperature  gradient  in  liquid 

k  a  constant 

Kg  effective  distribution  coefficient 

K  equilibrium  distribution  coefficient 
o 

L  moles  of  liquid 

L  initial  moles  of  liquid 

m  slope  of  liquidus  line 

n  dimensionless  distance  from  interface  into  liquid 

q.  fraction  frozen  of  pure  component 


8  moles  of  solute  remaining  in  the  ^stem  after  a  certain  degree  o£ 

solidification 

Sq  total  noles  of  solute  in  the  s7stem 

Tg  equilibrium  temperature 

T  initial  temperature 

W  weight  fraction  of  solute  in  euteotio 

V  weight  fraction  in  liquid  a  distance  x  from  interface 

Wj  mean  bulk  weight  fraction  in  liquid 

W  weight  fraction  of  solid 

8 

W  initial  weight  fraction  in  liquid 
o 

X  distance  from  initial  point  of  freezing 

x'  distance  from  interface 

x^  initial  mole  fl-action  of  liquid  phase 

X  mole  fraction  of  liquid  phase 

a  Kf/D 

$  boundazy  layer  thickness 

<  dimensionless  eatectic  composition 


INTRODUCTION 

Normal  freezing  is  a  term  applied  to  the  progressive  freezing  of  a 
liquid  charge.  It  is  veil  known  that  -atiea  a  solution  undergoes  a  normal 
fl-eezing  process  a  nonlinear  composition  profile  will  result. 

The  first  quantitive  explanation  of  this  phenomena  vas  developed  by 
Pfann  (1).  Historically,  normal  freezing  has  been  applied  to  the  purifi- 
cation of  materials  exhibiting  solid  solution  phase  behavior;  therefore, 
Pfann  approximated  the  solid  solution  phase  diagram  with  straight  lines  as 
i^oun  in  Figure  1.  From  thermodynamics  it  can  be  shown  that  this  approximation 
becomes  better  as  the  solute  concentration  decreases.  The  assumption  of  lovr 
solute  concentration  is  not  a  serious  restriction  as  usually  normal  freezing 
is  applied  to  Just  such  cases.  It  is  convenient  to  define  the  equilibrium 
distribution  coefficient  as  shown  below. 

This  coefficient  can  be  considered  a  constant  because  of  the  straight  line 
assumption  made  for  the  p^ase  diagram.  Other  assumptions  made  by  Pfann  are: 

1.  Diffusion  of  solute  in  the  solid  is  negligible 

2.  Mixing  in  the  liquid  phase  is  complete* 
By  definition 

Cg  »  Va  (1) 

The  liquid  phase  composition  is  given  by 

Ca.«  s/(l-g)  (2) 

Combining  Equations  (1)  and  (2) 

Cg  =  K^s/d-g)  (3) 

assuming  fireezing  to  take  place  in  a  tube  of  unit  volume.  Let  an  element  of 


volume  dg  fi-eeze.  The  concentration  of  this  frozen  volume  is: 

C    =  -ds/dg  C^) 
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Eliminating  Cg  from  Equation  (3)  and  expressing  as  an  integral: 
/  f  =  -  Ko  ^*  dg/(l-g) 


(5) 


o 

Integrating: 


s=So(l-g)^^  (6) 

ELirainating  s  by  using  Bquation  (3)  and  noting  that  C^  =  s^ 

This  is  Pfann's  expression  for  solute  distribution  in  a  solid  ^en  the  liquid 
phase  is  completely  mixed. 

Tiller,  Jackson,  Rutter,  and  Chalmers  (2)  have  derived  the  expression 
for  solute  distribution  assuming  that  liquid  phase  mixing  is  negligible  and 
mass  transfer  is  caused  only  by  diffusion.  As  in  Pfann's  analysis,  solid 
phase  diffusion  is  assumed  to  be  negligible  and  Kq  is  assumed  to  be  constant. 
For  sake  of  analysis  the  fl-eezing  is  assumed  to  take  place  uniaxiail.ly  in  a 
tube  of  unit  cross-section.  As  freezing  occurs  solute  is  rejected  by  the 
solid  and  diffuses  into  the  liquid  causing  a  build  up  os  solute  in  the  liquid 
near  the  interface  and  hence  a  concentration  gradient.  A  sketch  of  this 
concentration  gradient  is  presented  in  Figure  2.  The  rate  of  freezing  will 
be  constant  and  equal  to  f.  Assuming  steady  state,  the  flows  into  and  out 
of  the  liquid  volume  element  dx  are; 
Flow  into  element 


Bulk  Flow         f  [c  +  (^,)  dx'J 


Diffusion         "  ^  S» 


V 
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Flow  out  of  element 

Bulk  flow 

to 

DifltLsion 

-=Li^<«  +  €.>^'>] 

Assuming  no  accumulation 

Flow  in  -  flow  out  B  0 

fL''*(i.J'^'J-''iT- 

•^o  +  4i-*sr€.5^ 

•J   (8) 

Simplification  gives: 

' 

..     ^i--5"° 

(9) 

The  boundary  conditicms  are: 

ac'    •    Cl  «  Co 

X»  a  0     C-  =  C  /K 
L    0  0 

The  solution  is  then: 

fic» 

Ct  =  A  e    ^  +  C^ 
L             0 

(10) 

For  steady  state  it  is  necessaiy  that  C  =  C^,  also  it  is  apparent  that: 

A  =  C^  /  K„  -  C« 
0  '   0    o 

(11) 

combining  Equations  (10)  and  (11): 

r   r   1  1     ®     ^ 

• 

^L  *  ^0  L^  +   K   ® 

0 

This  is  the  composition  profile  as  it  exists  in  the  liquid  phase. 

Since 

"we  desire  the  solid  phase  profile  oui 

'  analysis  is  not  c<xaplete. 

Figure  3 

shows  a  sketch  of  what  the  distribution  profile  should  look  like  if  it 

follo\TOd  the  expected  general  conditions. 

1.  The  initial  solidification  would  have  a  value  of  K^C^ 

2.  C  should  be  approached  asymptotically  as  the  liquid 
phase  disappears. 

3.  C  must  rise  continuously  firom  K^C^jto  Cq  tram,  a  solute 
material  balance  it  is  seen,  Figure  3t  that  the  area  between  C^  and  Cg 
curves  must  equal  the  area  betxreen  the  Ct  and  C^  curves.  This  can  be  re- 
presented analytically  as  Equation  (13)* 

S^  (C^  -  Cg)  dx  =  f     (C^  -  C^)  dx'  (13) 

If  it  is  assumed  that  the  rate  of  approach  of  C.  to  Cq  with  distance  is 
proportional  to  (C  -  C  )»  then: 

cl(Co  -  C-) 

—V-^  «  -  or(C^  -  C3)  (lif) 

Equation  (14)  satisfies  the  above  set  of  conditions.  The  solution  to  this 
equation  is: 

C  -  C  =  A  e*^  +  B  (15) 

OS  '  -" 

applying  the  previously  mentioned  conditions: 

(Co   -   Cg)  0  X  • 

s   o  o 
the  solution  is: 

Cfl  -  Co  [1  -  (1  -  Ko)  e-  *^]  (16) 

substituting  Equations (12)  and  (16)  into  (I3)  gives: 

fjc' 

°o^^  "  V  f  ^"  "*  <^  '  °o  "4^  .^*  «"  "^  ^'         ^17) 
.0  o   0 

IMs  equation  yields 


a3-|-  (18) 

using  this  value  of  a,  Biuation  (l6)  becomes: 

Cs  =  ^o  t  -  (l-^o)  ^   (-  T  ^)J  ^^9^ 

The  derivation  of  this  equation  makes  no  provision  for  the  excess  of  solute 
that  will  appear  in  the  tenninal  region  and  therefore,  cannot  be  applied  to 
the  final  stages  of  solidification.  It  has  been  shown  by  Smith,  Tiller  and 
Ratter  (3)  that  the  tenninal  transient  will  appear  as  shown  in  Figure  4. 
The  equation  for  the  tenninal  transient  is  quite  cumbersome  and  idll  not  be 
included  here* 

Bie  first  two  equations  (7)  and  (19)  represent  extremes  in  the  assumed 
liquid  phase  mixing.  The  analysis  by  Burton,  Prim  and  Slichter  (4),  represents 
the  condition  of  partial  mixing  by  using  a  boundary  layer  approach.  Using 
Equation  (9)  and  applying  the  boundaiy  condition: 

C  «  C^j    at    x>  »  6 
dis  assumes  that  the  concentration  gradient  exists  only  in  a  boundary  layer 
of  thickness  6  and  that  outside  the  boundary  layer  the  concentration  is  Ct  . 
From  the  steady  state  analysis  it  is  seen  that  another  boundary  condition  is: 

(C^  -  Cg)f  +  D  ^t  =  0    at    x«  a  0  (20) 

applying  these  boundary  conditions  to  Equation  (9)  gives: 
C-  -  C-   -f6/D 

Defining  a  new  variable  KL   :    ,.. 

called  the  effective  distribution  coefficient  and  substituting  K  and  K« 
into  Equation  (21)  gives  in  logarithmic  fona. 
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!Ihis  is  the  working  form,  of  the  equation  used  in  this  investigation.  Notice 
this  eqvtation  is  consistent  vith  the  extreme  cases  of  no  mixing  and  complete 
mixing: 

No  mixing  -«—«-—  few  and  Kg  =  1 
complete  mixing  — —  f  «  0  and  Kg  =  K^ 

Ihe  qualitative  representation  of  solute  distribution  along  the  growth 
axis  is  represented  in  Figure  5t  for  the  three  cases  discussed  here.  It  is 
seen  that  Pfann's  complete  mixing  model  gives  the  best  segregation  xAile  the 
no-mixing  model  of  Tiller  et.  al.  gives  the  poorest  segregation.  Therefore, 
mixing  of  the  liquid  phase  would  be  expected  to  increase  segregation. 

All  the  previous  discussion  has  concerned  mixtures  having  solid  solution 
phase  behavior.  However,  this  thesis  is  concerned  with  binaiy  organic  systeus 
t&ich  in  general  exhibit  simple  eutectic  phase  behavior.  A  typical  phase 
diagram  for  such  an  oi^anic  system  is  shown  in  Figure  6.  Notice  that  the 
solid  phase  obtained  from  an  organic  mixture  of  composition  C_  would  be  ex- 
pected to  consist  of  pure  component  A,  in  other  words  the  value  of  Kq  is  jsero. 
This  means  that  technically  the  previous  equations  applying  to  solid  solutions 
are  not  applicable  to  eutectics  as  a  review  of  the  derivations  will  show. 

For  complete  mixing  of  a  liquid  phase,  simple  eutectic-foraing  mixture 
the  analysis  is  easy.  Pure  solute  will  freeze  out  until  the  liquid  phase 
reaches  the  eutectic  composition;  then  the  solid  phase  will  freeze  at  the 
eutectic  composition,  see  Figixre  ?• 

Wilcox  (5)  has  analysed  and  obtained  an  equation  for  simple  eutectic- 
forming  mixtvires  assuming  a  boundary  layer  exists  next  to  the  interface  in 
■tdiich  no  mixing  occurs,  the  bulk  of  the  liquid  is  assumed  well  mixed. 


According  to  micox,  a  certain  fraction  of  liquid  will  freeze  producing  pure 
component  in  the  solid  i^ase.  The  expression  for  this  fraction  is: 

(ee-^  -  1)  /ee-^  (23) 
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Vhen  more  than  this  fraction  of  liquid  is  solidified  the  eutectic  composition 
is  reached  in  the  boundary  layer  and  the  distribution  is  given  by  Equation  (2k) 

W^  =  W^  {l -[§E|^Jaxp  (e-^/1  -  e-^)}  (2^) 

Figure  8  shows  the  calculated  solute  profile  obtained  by  applying  Wilcox's 
equations  to  a  hexadecane-benzene  system  containing  3  mole  i  hexadecane. 

Past  experimental  work  with  solid  solutions  has  shown  the  actual 
segregation  to  be  less  than  that  predicted  by  theoiy.  The  chief  reason  for 
this  lack  of  verification  is  thought  to  be  due  to  a  phenomenon  first  proposed 
by  Rutter  and  Chalmers  (6)  termed  "constitutional  subcooling".  As  the  inter- 
face of  a  freezing  solution  rejects  solute  a  concentration  gradient  ^dll 
exist  Joi  the  liquid  as  already  shown  in  Figure  2.  Since  the  equilibrium 
freezing  temperature  depends  upon  concentration  (as  determined  by  the  phase 
diagram),  an  equilibrium  temperatiire  gradient  will  exist  similar  to  curve  "a" 
of  Figure  9.  The  actual  temperature  gradient  imposed  in  the  liquid  can  be 
represented  by  curve  "b".  At  every  point  \&ere  curves  "b"  is  below  curve  "a" 
the  liquid  will  be  subcooled.  5y  representing  curves  "a"  and  "b"  analytically. 
Tiller,  et  al.  (2)  were  able  to  show  that  the  critical  rate  above  vixich 
subcooling  will  exist  is  given  by: 

£„.—-^ (25) 

(^)  (l-K^) 
^o 

This  equation  implies  that  there  is  a  minimum  rate  of  freezing  below  liiich  no 
subcooling  will  exist.  This  then  is  the  phenomena  knox«i  as  constitutional 
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supercooling  and  is  responsible  for  the  irreularities  observed  at  the  inter- 
face of  freezing  mixtiires* 


n 


PROCEDURE 

Die  apparatus  used  for  conducting  normal  freezing  experiments  consisted 
of  a  refrigerated  methyl  alcohol  bath  and  a  drive  mechanism  to  control 
fteezing  rates.  Ihe  bath  teii5>erature  was  varied  from  -10  C  to  -33  C  for 
different  runs,  and  was  controlled  so  that  the  maximum  variation  for  any 
teaperatixre  setting  was  approximately  one  degree.  To  keep  the  bath  at  a 
uniform  temperature  it  was  agitated  by  a  l/l2  H.  P.  Fultork  Labmotor  operating 
a  stirrer. 

The  drive  mechanism  consisted  of  a  long  upright  shaft  threaded  to  screw 
into  a  threaded  gear.  Keeping  the  gear  stationary  in  the  horizontal  plane 
and  allowing  the  gear  to  rotate  caused  the  shaft  to  ascend  or  descend  depend- 
ing on  the  direction  of  rotation.  A  rachet  assembly  powered  by  a  1/30  H.  P. 
motor  caused  the  gear  to  rotate.  The  rod  could  descend  at  a  constant  rate 
between  0.1  in./hr.  and  3«0  in./hr.  depending  on  the  rachet  setting.  The 
drive  was  equipped  with  an  automatic  shut  off  to  stop  the  motor  at  the 
completion  of  each  run.  To  allow  the  motor  to  run  might  damage  the  drive 
mechanism. 

In  making  a  typical  expeidmental  run  a  glass  tube  sealed  at  one  end 
was  filled  with  approximately  50  grams  of  the  specimen  solution.  The  vreight 
of  the  specimen  was  recorded  as  a  mass  balance  check.  The  tube  was  held  in 
place  vertically  over  and  just  touching  the  bath  by  means  of  a  clamp  attached 
to  the  descending  rod.  The  motor  was  started  and  the  tube  slowly  descended 
into  the  bath.  Usually  it  would  take  several  minutes  for  the  liquid  to 
start  freezing.  If  fl^ezing  had  not  occurred  in  the  first  several  minutes 
the  tube  was  tapped  gently  so  as  to  initiate  the  freezing.  After  the  entire 
specimen  had  frozen  the  drive  was  stopped. 
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nil  order  to  determine  the  composition  profile  it  was  necessary  to 
divide  the  solid  specimen  into  a  number  of  portions,  (usually  eight). 
Removal  of  each  portion  was  accomplished  by  raising  the  tube  out  of  the 
bath  so  the  volume  of  solid  to  be  removed  was  above  the  level  of  the  bath, 
this  portion  was  then  allowed  to  melt.  Generally  a  half  hour  was  required 
to  complete  the  melting  of  each  portion.  Care  was  exercised  at  this  point 
as  the  interface  tends  not  to  be  flat.  Upon  standing  a  time,  (depending 
on  tube  diameter)  the  interface  became  flat.  Only  a  melt  with  a  flat  inter- 
face will  give  an  accurate  composition  measurement.  Each  melted  portion  vras 
transferred  to  a  weighing  bottle  and  weighed  to  a  l/lO  milligram. 

A  Bausch  and  Lomb  reffactometer  was  used  to  measure  the  refractive 
index  of  each  increment.  A  composition  curve  versus  refractive  index  had 
been  determined  frcaa  mixtures  of  known  composition  and  xras  used  to  convert 
refractive  indexes  to  compositions.  All  the  refractive  indexes  were  read 
at  25.0  °C  to  an  acciiracy  of  +  0,0001.  The  calibration  curve  for  the 
hexadecane-benzene  system  is  given  in  appendix  (1). 

Three  liquid  systems  were  investigated,  they  were: 

benzene-normal  heptane 

benzene-cyclohexane 

benzene-hexadecane        ' 
All  starting  mixtures  of  benzene-hexadecane  contained  approximately  97  mole  ^ 
benzene.  The  heptane  gystem  contained  98  mole  i  benzene  Tjhile  the  cyclo- 
hexane  gystem  contained  96  mole  ^  benzene. 

For  some  runs  it  was  necessary  to  remove  dissolved  gas  from  the  specimen. 
This  was  accomplished  by  submerging  the  entire  tube  in  the  bath  until  the 
contents  has  became  solid,  then  evacuating  the  tube  and  allowing  the  solid 
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to  melt,  VJien  this  process  was  repeated  three  times  the  specimen  was 
effectively  degassed  and  no  bubble  formation  during  normal  freezing  was 
observed.  To  be  certain  the  degassing  process  did  not  redistribute  the 
solute  the  liquid  was  carefully  agitated  before  each  run. 

Mixing  of  the  liquid  portion  of  the  sample  during  freezing  was  accom- 
plished by  attaching  a  flat  disc  to  a  line  and  suspending  the  disc  horizontally 
in  the  liquid  phase.  Allowing  the  disc  to  reciprocate  up  and  down,  by  means 
of  a  motor,  caused  mixing  in  the  liqviid  phase.  Diis  arrangement  was  used 
in  favor  of  a  blade  stirrer  as  the  amovmt  of  agitation  caused  by  a  stirrer 
is  dependent  upon  the  amount  of  liquid,  (constantly  changing  during  a  run). 
Bie  mixing  caused  by  the  chosen  arrangement  was  considered  to  be  independent 
of  the  quantity  of  liquid. 

To  obtain  temperature  distributions  during  freezing  a  thermocouple  was 
suspended  in  the  liquid  at  a  fixed  distance  relative  to  the  tube.  As  the 
solid-liquid  interface  moved  up  the  tube,  temperature  measuraaents  were  taken 
at  various  times  so  that  a  temperature  profile  was  obtained.  Ihe  point  at 
viiich  the  interface  touched  the  thermocouple  was  recorded.  Temperatures 
were  read  to  an  accuracy  of  0.05  °C. 

In  order  to  better  vmderstand  the  solute  rejection  mechanism  photographic 
studies  were  made  of  the  solid-liquid  interface.  A  microscope  equipped  with 
a  camera  was  mounted  a  few  degrees  above  the  horizonal  position.  With  this 
arrangement  a  profile  microphotograph  of  the  interface  could  be  obtained.  In 
order  to  prevent  distortion  the  interface  was  photographed  in  a  tube  with 
polished  flat  slides.  All  photographs  were  made  at  the  same  magnification  of 
14  times.  Exposure  time  was  15  seconds;  this  long  exposure  time  required 
that  the  drive  mechanism  be  stopped  just  prior  to  exposure  to  insure  against 
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movement  of  the  tiibe.  A  fine  grain  film  (Kodak  Contrast  Process  Ortho) 
was  used  for  all  liiotographs  and  developed  ijith  a  high  contrast  developer 
(Kodak  D-11). 
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RESULTS 

Die  effect  of  the  followiiig  parameters  on  segregation  by  normal  freezing 
iaas  investigated. 

1.  Bath  temperature. 

2.  Mxing  of  liquid  please. 

3.  Degassed  specimen. 
k.       Tube  diameter. 

It  >as  found  that  each  individual  run  vxas  well  represented  by  the 
incomplete  mixing  model.  Figures  10,  11,  12  and  13  are  typical  of  the 
composition  profiles  obtained.  Wilcox's  equations  for  eutectic  systems  were 
found  not  to  apply  as  at  no  time  was  100  $2  pure  material  obtained.  Equation 
(7)  for  the  incomplete  mixing  model  can  be  put  in  the  follotdng  form. 

In  x/xq  =  (1  -  Kg)  In  Lq/L  (26) 

It  was  found  that  a  log-log  plot  of  x/x^  versus  Lq/L  for  all  runs  was  re- 
presented \iell  with  a  straight  line.  Figures  14-  and  15  are  typical  of  this 
correlation.  Because  of  the  excellent  agreement  it  was  possible  to  calculate 
Kg  from  the  slopes  of  these  plots.  Table  1  contains  values  of  the  distribu- 
tion coefficient  for  all  runs.  In  order  to  correlate  K«  and  f ,  Burton, 
Prim,  and  Slichter's  boundaiy  layer  model  was  employed.  From  Eqixation  (22) 
it  is  seen  that  a  semi-log  plot  of  (l/Kg  -  1)  versus  f  would  be  expected  to 
produce  a  straight  line  of  slope  -  6/d  and  intercept  of  (I/k^  -  1).  Previously 
it  was  found,  (?)  that  this  correlation  applied  well  to  the  cyclohexane- 
benzene  and  heptane-benzene  systems  as  shown  in  Figures  16,  17,  and  18.  The 
hexadecane-benzene  system.  Figures  19 t  20  and  21  xfere  best  represented  by 
two  straight  lines.  However,  a  single  straight  line  was  found  to  fit  the 
stirred  and  degassed  runs,  Figures  22  and  23. 
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Effect  of  Parameters 

Bath  teiaperattire 

The  effect  of  bath  temperature  on  nonstirred  runs  can  be  seen  from 
Figure  20.  Five  data  points  were  taken  for  a  low  -30  °C  bath  ten^Jeratxire 
and  seven  points  were  taken  at  a  higher  -18  °C  bath  temperature.  It  is  seen 
that  over  this  temperature  range,  bath  temperature  appears  to  have  no  effect 
on  segregation.  Figure  22  shows  the  effect  of  bath  temperature  on  a  stirred 
system;  it  is  seen  that  segregation  for  the  hexadecane-benzene  system  decreases 
xdth  decreasing  bath  temperature. 

I-Iixing  of  liquid  phase 

For  the  cyclohexane-benzene  system  it  was  found  that  mechanical  stirring 
of  the  liquid  phase  increased  segregation  as  illustrated  in  Figure  2k.     Hotrever, 
Figure  25  shows  that  similar  stirring  of  the  hexadecane-benzene  system  produced 
less  segregation  than  was  obtained  for  unstirred  runs. 

Degassed  specimen 

All  runs  produced  gas  bubbles  at  the  interface  except  \^en  the  specimen 
was  purposely  degassed  as  mentioned  in  the  procedure.  The  effect  of  gas 
bubbles  was  investigated  only  for  the  hexadecane-benzene  system.  Figure  26 
shows  that  removal  of  the  dissolved  gas  decreased  segregation. 

Tube  diameter 

The  effect  of  tube  diameter  was  investigated  for  the  hexadecane-benzene 
system.  Figure  2?  is  a  summary  of  the  results  that  were  presented  in  Figures 
19.  20,  and  21.  It  is  seen  that  segregation  increases  with  increasing  tube 
diameter. 
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Photographic  Studies 

Micropihotographs  were  taken  of  the  interface  of  pure  benzene  and  the 
hexadecane-benzene  systea.  These  photographs  are  presented  as  plates  (1) 
through  (5).  It  is  observed  that,  at  freezing  rates  betvieen  1.0  and  2,5  in./hr. 
the  Interface  maintains  an  irregular  dendritic  surface.  At  freezing  rates 
less  than  1.0  in./hr.  the  iixterface  shows  a  change  from  dendritic  formation 
to  a  more  regular  but  corrugated  surface. 

It  was  difficult  to  photograph  the  interface  of  pure  benzene  as  this 

substance  is  not  opaque  \Aen  frozen,  as  is  the  hexadecane-benzene  system,  but 

remains  clear.  Plate  (6)  is  the  best  iSiotograph  that  could  be  obtained  for 

pure  benzene.  Die  interface  is  faint  but  discernible  as  a  sharp  line  across 

these  photographs.  The  white  blotches  in  the  background  above  the  interface 

are  reflections  from  bubbles  still  attached  to  the  interface. 

■  7   r 

Constitutional  Subcooling 

Temperature  profiles  xrere  taken  in  order  to  solve  Equation  (25)  and  to 

determine  if  constitutional  subcooling  was  occuring.  A  typical  profile  is 

presented  in  appendix  (2).  !Ihe  rate  of  freezing  was  found  to  have  a  small 

effect  on  the  temperature  gradient.  Equation  (25)  revealed  that  all  runs 

were  subject  to  constitutional  subcooling  effects. 
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DISCUSSION 

Bath  temperature 

It  might  be  expected  that  bath  temperatxire  would  not  effect  segregation 
as  the  position  of  the  interface  relative  to  the  bath  level  idll  vaiy  so  that 
the  interface  is  at  the  freezing  temperature  and  not  necessarily  at  the  bath 
temperature.  This  was  found  true  for  \mstirred  runs.  The  effect  of  bath 
temperature  on  stirred  runs  is  thought  to  be  due  to  the  rather  high  eutectic 
temperature,  (see  appendix  3)  of  the  hexadecane-benzene  system.  It  was 
observed  that  for  the  same  bath  temperature  the  interface  for  stirred  runs 
was  significantly  lower  telative  to  the  bath  than  for  unstirred  runs.  This 
suggests  that  the  interface  is  not  at  the  equilibrium  temperature,  but  is 
below  this  temperature.  For  the  hexadecane-bezene  system  if  the  interface 
temperature  were  only  a  few  degrees  below  the  equilibrium  temperature  it 
would  be  possible  to  freeze  liquid  at  the  eutectic  composition  resulting  in 
lowered  segregation. 

Mixing  of  liquid  phase 

laxing  of  the  liqviid  phase  for  the  cyclohexane-benzene  system  improved 
segregation  as  expected,  since  the  eutectic  teu^ierature  for  this  ^stea  is 
well  below  the  bath  temperature.  The  reason  for  decreased  segregation  of 
the  stirred  hexadecane-benzene  system  is  believed  to  be  due  to  the  high 
eutectic  point  of  this  system  as  mentioned  above. 

Degassed  specimen 

The  decrease  in  segregation  caused  by  degassing  of  the  specijnen  is  in 
keeping  with  the  boundary  layer  theory.  It  would  be  expected  that  rising 
bubbles  of  gas  would  cause  a  certain  degree  of  mixing  to  occur  in  the  liquid. 
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Sine©  mixing  increases  segregation,  the  removal  of  the  source  of  this  mijcing 
should  cause  a  drop  in  segregation.  It  should  be  pointed  out  that  this 
conclusion  is  not  necessarily  inconsistent  with  the  lack  of  separation  ob- 
served in  the  stirred  system  as  the  agitation  caused  by  rising  gas  bubbles 
is  very  mild  compared  to  the  vigorous  forced  convection  caused  by  the  stirrer. 
The  interface  distance  was  not  observed  to  change  because  of  degassing, 
indicating  the  effect  of  nonequilibrium  freezing  was  not  present. 

Tube  diameter 

Perhaps  the  strongest  support  for  the  boundary  layer  model  is  found  in 
the  effect  of  tube  diameter  on  segregation.  From  elementary  boundary  layer 
theory  it  is  known  that  as  convection  in  a  fluid  increases  the  stagnant 
boundary  layer  decreases.  Also  increasing  tube  diameter  would  be  expected 
to  result  in  an  increase  of  natural  convection.  As  shown  in  the  results 
the  largest  diameter  tube  gave  the  better  segregation  at  a  particular  freezing 
rate  and  the  smallest  tube  gave  the  poorest  segregation.  Furthermore,  it  was 
noticed  that  for  rates  less  than  1.0  in./hr.  the  smallest  tube  size  had  the 
greatest  slope,  Figure  2?  lAile  the  largest  tube  had  the  least  slope.  The 
slope  of  these  lines  as  given  by  Equation  (22)  is  -  6/D.  If  it  is  assumed 
that  D  remains  unchanged  then  6  must  decrease  as  the  tube  size  increases, 
thus  satisfying  the  above  mentioned  expectations.  The  branch  of  the  hexadecane 
line.  Figure  2?  for  rates  greater  than  1.0  in./hr.  does  not  satisfy  the 
requirement  that  the  smallest  tube  should  give  the  greatest  slope,  this  may 
be  explained  as  follows. 

Constitutional  subcooling 

The  interface  of  the  hexadecane-benzene  system  was  foxand  to  possess  an 
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irregular  surface.  However,  the  microphotograjJis  reveal  a  change  in  the 
type  of  surface  irregularities  occurring  at  a  rate  of  1.0  in./hr..  For 
rates  is  excess  of  1.0  in./hr.  the  interface  becomes  qviite  dendritic,  see 
Plates  W   and  (5).  It  is  believed  that  the  discontinuous  slope  for  the 
hexadecane-benzene  curve  shown  in  Figure  2?  is  a  direct  consequence  of  this 
change  of  interfacial  morphology.  Both  changes  appear  to  occur  at  a  freezing 
rate  of  1.0  in./hr..  The  nature  of  the  dendritic  interface  is  such  that  it 
would  appear  that  trapping  of  liquid  could  readily  occur  for  freezing  rates 
greater  than  1.0  in./hr..  Any  trapping  of  liquid  would  be  expected  to  lessen 
the  degree  of  segregation  and  hence  a  change  in  slope  of  the  line  in  Figure 
2k,    For  freezing  rates  less  than  1.0  in./hr.  the  interface  took  on  a  more 
corrugated  appearence.  Plates  (1),  (2)  and  (3).  The  amount  of  liquid  trapping 
for  this  type  of  interface  would  be  expected  to  be  much  less  than  for  the 
dendritic  interface.  Nevertheless  some  degree  of  liquid  trapping  could  be 
expected,  thus  accounting  for  the  non-zero  equilibrium  distribution  coeffi- 
cients obtained  as  the  exprapolated  intersection  of  the  plot  of  Figure  21+ , 
at  zero  freezing  rate.  !Ihis  transition  from  a  corrugated  to  a  dendritic 
interface  is  believed  to  be  analogous  to  the  transition  from  a  celluar  inter- 
face to  a  dendritic  interface  as  observed  in  metal  systeas  by  other  investi- 
gators (8),  (9),  and  (10). 

Furthermore,  the  better  segregation  obtained  for  the  cyclohexane-benzene 
and  the  heptane-benzene  systems  is  believed  to  be  due  to  a  celluar  rather 
than  a  dendritic  interface  for  these  systems  over  the  freezing  rates  investi- 
gated. It  can  be  shown  that  the  transition  from  a  celluar  to  a  dendritic 
interface  would  be  expected  at  a  lower  freezing  rate  for,  the  hexadecane- 
benzene  system.  Both  the  solid  solution  and  eutectic  boundary  layer  models 
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account  for  this  fact.  For  solid  solutions  Equation  (25)  gives  the  critical 
value  of  f  above  vhlch.   subcooling  occurs.  This  can  be  rearranged  to: 

o 

Using  Wilcox's  equations  for  eutectic  systems  the  critical  f  is  given  by, 

appendix  (^): 

G  ,  (i)  e  fiCf/D)  (28) 

Xn.  factors  in  Equations  (26)  and'  (27)  are  constant  except  D  and  f ;  However, 
the  ratio  D/f  must  also  be  constant,  therefore: 

D  =  kf  (29) 

Strictly  speaking  Equations  (27)  and  (28)  apply  only  at  the  transition  from 
a  flat  interface  to  a  celluar  interface;  however,  similar  considerations 
would  be  expected  to  be  valid  at  the  transition  of  a  celluar  to  a  dendritic 
interface.  The  table  below  shows  the  values  of  D,  (U). 

D  (cm  /sec)     6  (cm) 
^clohexane-benzene        1.3^  x  10"^   4.2  x  10"^ 
heptane-benzene  1.13  x  lO"'^   4.2  x  10"^ 

hexadecane-benzene         0.70  x  10~-^   4.7  x  10~^ 
Because  of  the  lower  diffusion  coefficient  the  transition  would  be 
expected  to  occur  at  a  lower  freezing  rate  for  the  hexadecane-benzene  system. 
As  assumed  above  and  as  expected  from  boundaiy  layer  theozy  the  boundary 
layer  thickness  should  be  independent  of  the  natxire  of  the  solute  and  of  the 
freezing  rate.  The  above  table  shows  6  as  calculated  from  the  slope  of  the 
curves  in  Figures  16,  18,  and  20  and  from  the  value  of  D.  The  boundary  layer 
thickness  appears  to  be  relatively  constant,  thus  adding  support  to  the 
boundary  layer  theory. 


22 


CONCLUSIOKS 

G3ie  variation  of  segregation  by  normal  freezing  vdth  freezing  rate 
is  well  represented  by  Burton,  Prim  and  Sliohter's  incomplete  mixing  boundary 
layer  model.  The  effects  of  mixing,  tube  diameter  and  bath  temperature  can 
be  explained  in  terms  of  this  boundaiy  layer  model.  The  failure  of  Tdcox's 
model  for  eutectic  systems  is  thought  to  be  due  to  constitutional  subcooling 
which  in  t\im  results  in  the  trapping  of  liqviid  in  the  solid. 
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Run  No, 


Table  1 

Freezing  Rate 
(in./hr.) 

Distribution 
Coefficient 

Bath 
Temperature 

Hexadecane-Benzene 

3A"  Dia. 

Tube 

2.25 

.733 

-  20 

1.75 

.721 

-"20 

1.25 

.668 

-  20 

.75 

.6^ 

-  20 

3.0 

.807 

-  20 

.50 

.623 

-  20 

1.75 

.717 

-33 

2.75 

.789 

-33 

.50 

.^37 

-  33 

.75 

.627 

-33 

.50 

.607 

-33 

.25 

.525 

-  20 

•10 

.485 

-  20 

1 

l^ 

5 
6 

7 
8 

9 

U 
12 

13 
14 
68 
71 

Hescadecane-Benzene     1.0"  Dia.  Tube 

22  1.75  .684  -  20 

23  2.75  .716  -  20 

36  2.5  .709  -  20 

37  1.75  .694  -  20 

38  1.25  .635  -  20 

39  .75  .635  -  20 
42           3.0             .714  -  20 
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Table 

1  cont'd 

Run  No. 

Freezing  Rate 
(in./hr.) 

Distribution 
Coefficient 

Bath 
Temperature 

^3 

1.0 

.636 

-  20 

69 

.25 

.538 

-  20 

Hescadecaae-Benzene 

1/2"  Dia. 

Tube 

57 

2.75 

.771 

-  20 

58 

2.75 

.785 

-  20 

59 

2.0 

.747 

-  20 

60 

1.0 

- 

.689 

-  20 

^ 

1.5 

.760 

-  20 

62 

.50 

.593 

-  20 

63 

2.75 

.799 

-  20 

6k 

.75 

.641 

-  20 

65 

1.25 

.717 

-  20 

66 

1.00 

.676 

-  20 

67 

.25 

.520 

V-  20 

Degassed     Heocadecane-Benzene 

3/4"  Dia. 

Tube 

kl 

3.0 

.758 

-  20 

40 

1.75 

.756 

-  20 

45 

1.75 

.796 

-  20 

46 

1.75 

.829 

-  20 

47 

1.75 

.795 

-  20 

46 
49 

2.75 
2.75 

.859 
•860 

-  20 

-  20 
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Table 

1  cont'd 

Run  No. 

Freezing  1 
(In./hr 

Rate 

Distribution 
Coefficient 

Bath 
Temperature 

50 

.75 

.760 

-  20 

51 

1.25 

.711 

-  20 

52 

1.25 

.802 

-  20 

53 

1.0 

.810 

«  20 

^    '. 

.50 

.753 

-  20 

,  55 

.50 

.708 

-  20 

56 

.50 

.715 

-  20 

Stirred  Hexadecane-Benzene     3A"  Dia.  Tube 

17  2.75  .809  -  18 

18  .50  .681  -  18 

19  1.75  •757  -18 

20  .75  .698  -  18 
a  2.25  .790  -  18 

15  1.25  .8105  -  28 

16  2^00  .8125  -  28 
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SOLUTE  CONCENTRATION 


Fig^.  1,         Portion  of  solid  solution  phase  diagram. 
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Fig  ,  2.,     Distribution  of  solute  ahead  of  interface, 
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Fig  .  3 ,  Bistribution  of  solute  in  solid  and  liquid  phases. 
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Fig  •  '^.   Concentration  profile  of  frozen  specimen . 
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Fig  •  5,  '^-Solute  distribution  along  growth  axis, 
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Fig  .  6.         Tjrpical  simple  eutectic  phase  diagram. 
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Fig  ,  7.  Solute  distribution  profile  of  simple  eutectic 
forming  mixture. 
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Fig  .  8.  Calculated  solute  distribution  using  Equation  (22) 
for  3  mole  %  hexadecane  and  benzene.. 
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Fig  ,  0.  Constitutionaly  supercooled  zone  ahead  of  interface. 
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Fig,  10.     Composition  profile. 
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Fig  .  14.   LOG  xAo  versus  LGG  L  /L. 
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Fig.  15.  log;  X/Xq  versus    LOG  Lq/L, 


3U 


0.2. 


3/4  Inch  Tube  Diameter 
Cyclohexane-Benzene 


-4- 


0.5 


1.0  1.5  2.0 

FREEZING  RATE   (inchAour) 


2.5 


Fig.  15.      LOG((i/K£)  -1)  versus  freezing  rate. 
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Fig.  18,      LOGi((i/Kg)  .1)  versus  freezing  rate. 
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Fig.  19.       LOG((l/K£)  -1)  versus  freezing. 
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Fig.;  2(i    LOG  T(l/%')  -1)  versus  freezing  rate. 
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Fig.  21.      L0G((1/Ke)  -1)  versus  freezing  rate. 
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Fig  .  22  .   LOG  (Xl/Kg)  -1)  versus  freezing  rate. 
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Fig  .  23  .    LOG  ((1/Ke)  -l)  versus  freezing  rate. 
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Fig,  25,  -   L0G((1/Ke)  -1)  versus  freezing  rate. 
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PROFILE  OF  INTERFACE 

hexadecane-benzene 
rate:  0,$0  in./hr. 
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PROFILE  OF  INTERFACE 
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Appendix    ( l) 
Refractive  index  curve  for  hexadecane  in  benzene. 
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Appendix  (2) 
Temperature  profile. 


15.0 


o 


caio.O 


n 


5.0 


0.0, 


I 


Hexadecane-benzene 
Freezing  rate:   2.5  in/hr. 
Bath  temperature  -20  C 


I 


10 


15 


20         25         30         35 

nm  (KiisiUTES) 


40 


45 


Sh 


Appendix    (3) 

Phase  diagrams  of  hexadecane-benzene,  cyclohexane- 
benzene,  and  heptane-benzene. 
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Appendix  (4) 
VJilccK's  composition  profile  for  eutectic  systeas  (5)  is  givai  as: 


$        e 


Applying  the  following  conditions: 

$  s  w/w,    at  n 
o 

$#  =  W./wi   at  bulk  concentration 


Gives 

Substituting  for  n  and  b  gives 


•  (x-  -  ^A 


C  =  C.  e 

i 

Using  this  value  of  the  concentration  gives  the  equilibrium  temperature  as: 

Taking  the  derivative  with  respect  to  x  gives: 

6f 

(^)    =imC,e  ^    ■ 
dx'  x'=0   D   X 

For  constitutional  subcooling  to  exist: 

d3E'=^ 

After  rearranging  this  gives: 

6f 
0       f  .   D 

mC    *  D  ® 


NORML  FREEZING  OF  ORGANIC  LIQUIDS 


by 


DAVID  ARTHUR  IRVIN 
B.  S.,  San  Jose  State  College,  1961i 


AN  ABSTRACT  OF  A  MASTER'S  THESIS 


submitted  in  partial  fulfillmet  of  the 


requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Chemical  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 


1965 


Theoretical  correlations  predicting  the  segregation  obtained 
by  normal  freezing  of  binary  systems  have  been  derived  for  both 
solid  solutions  and  eutectic -forming  mixtures  by  other  investi- 
gators. It  has  been  observed  in  metal  systems  that  these  correla- 
tions do  not  accurately  predict  segregation.  It  is  believed  the 
poor  agreement  is  due  to  trapping  of  the  liquid  within  the  solid 
phase.  This  trapping  is  in  turn  believed  caused  by  a  phenomena 
called  constitutional  subcooling. 

The  composition  profiles  presented  in  this  thesis  were 
obtained  for  binary  organic  systems.  The  expression  derived  by 
Wilcox  to  apply  to  just  such  systems  (simple  eutectic)  failed 
to  predict  the  degree  of  segregation.  It  was  found  that  a  good 
correlation  was  obtained  with  Burton,  Prim  and  Slichter's  in- 
complete mixing  boundary  layer  model  even  though  it  was  derived 
for  solid  solution  type  systems.  The  effect  of  mixing  in  the 
liquid  phase,  tube  diameter,  and  bath  temperature  tend  to  support 
the  boundary  layer  model.  Furthermore,  it  was  noted  that  for  the 
three  organic  systems  investigated  the  calculated  thickness  for 
the  boundary  layer  appeared  not  to  vary  a  significant  amount  as 
required  by  boundary  layer  theory.  The  failure  of  the  equations 
applying  to  eutectic  systems  is  thought  due  to  constitutional 
subcooling.  It  was  found  that  the  entire  spectrum  of  freezing 


rates  investigated  were  subject  to  constitutional  subcooling 
effects.  From  photographic  studies  of  the  interface  two  types 
of  interfacial  geometries  were  identified,  a  dendritic  and  a 
corrigated  interface.  The  change  from  one  form  to  the  other 
appeared  to  occur  at  a  freezing  rate  of  1.0  in/hr,  with  the 
dendritic  form  occuring  at  rates  greater  than  1,0  in/hr.  The 
graphical  representation  of  the  data  showed  a  discontinuity  to 
occur  at  approximately  the  same  rate  as  the  transition  in  inter- 
facial morphology  (1,0  in/hr).  The  discontinuity  is  thought 
due  to  increased  liquid  trapping  caused  by  the  dendritic  nature 
of  the  interface. 


